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£^ . We report the first direct measurement of the proton 3-D Born-Oppenheimer (BO) potential in 

any material. The proton potential surfaces in the hydrogen bonded superprotonic conductor 
Rb3H(S04)2 are extracted from the momentum distribution measured using Deep Inelastic Neutron 
£ — , Scattering(DINS). The potential has a single minimum along the bond direction, which accounts for 

£NJ ■ the absence of the antiferroelectric transition seen in the deuterated material, and a saddle point off 

the bond direction for tunneling into the next well with a barrier height of 350 meV. The measured 
potential is in qualitative agreement with phenomenological double Morse potentials that have been 
■ used to describe hydrogen bonds in other systems. 
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The shape of the Born-Oppenheimer potential of the proton in hydrogen bonded systems plays an essential role 
in a variety of phase transitions, in proton conductivity in biologically and technically important materials, and 
in the mechanisms of a wide variety of chemically and biologically interesting reactions. As such, it has been the 
subject of extensive experimental, theoretical and phenomenological investigations over many years. It can sometimes 
be inferred from neutron and light scattering measurements of the excitations in the system, but these are rarely 
sufficient to provide an unambiguous characterization of the potential, particularly if the potential energy surface is 
highly anharmonic. Ab-initio calculations of the potential are subject to approximations of uncertain validity, and 
are usually assumed to be accurate when they agree amongst themselves. We present here the first direct, model 
independent, measurement of the of the 3-D proton Born-Oppenheimer potential in any system, and in particular, in 
■ a hydrogen bonded system of technological interest, Rb3H(S04)2- [? ] These measurements can be used to investigate 
the accuracy of ab-initio calculations and provide a data base for their improvement, to interpret neutron and light 
^*<0 ■ scattering measurements, and as input to simulations of the dynamics of the system containing the bond. In this case, 
the material is a superprotonic conductor, and the form of the potential is essential information in simulating the 
transport of the protons. The measurement agrees qualitatively with earlier phenomenological attempts to characterize 
the O-H-0 bond using empirical Morse potentials. [? ] A previous determination of the potential by means of inelastic 
neutron scattering and Raman scattering in the system we have measured[? ] is qualitatively wrong, indicating the 
difficulty of determining the potential from spectroscopy. The present measurements, made with neutron Compton 
scattering, have an rms error of about 20mev over most of the domain in which the potential can be measured. The 
technology exists to reduce that to 2mev. 

We extract the information needed to determine the potential from a measurement of the proton momentum 
distribution by a procedure suggested by Reiter and Silver. [? ] In order to do this, the proton must be at a center of 
inversion symmetry, the interaction with other protons must be negligible, and the temperature must be lower than 
any excitations of the proton vibrations. All of these conditions can be satisfied in RB3H(S04)2, which is of interest 
both as a superprotonic conductor at high temperatures and because of the large isotope effect it exhibits at low 
temperatures. The deuterated material has an anti-ferroelectric phase transition at 82K, while the protonated version 
has no anti- ferroelectric transition at all. [? ] An earlier experiment on KDP[? ] was able to extract the potential 
along the bond, but it was necessary to assume that the potential was separable because of the insufficient quality of 
the data. We make no assumptions here about the a priori form of the potential. 

The geometry of this H-bonded system consists of pairs of sulphate tetrahedra hydrogen bonded together, well 
separated from other pairs. The distance between nearest neighbor protons is 5.07A, leaving the protons on different 
bonds effectively isolated from each other. The A2/a space group places the average proton position at a center of 
inversion. This allows us to assume that the ground state momentum wavefunction is real. We will argue later that 
this is not just a statistical average, and that the hydrogen bond is itself symmetric. The potential can be expected 
to change with temperature, and indeed, we observe that it does, but we will only be concerned here with the low 
temperature region. 

Several experimental studies suggested the presence of a double well proton potential along the hydrogen bond [? 
? ? ] at same time that other experimental studies as well as theoretical calculations supported single well potential 
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models [???]. In Fillaux et al's work[? ], which investigated several related compounds, the peaks in INS spectra 
were assigned to a slightly asymmetric potential with a double well shape along the bond direction. We find that the 
potential is unambiguously a single well along the bond. There is a barrier to escape of the proton from the well into 
neighboring wells off the bond axis, with a barrier height of approximately 350 meV, in qualitative agreement with 
NMR experiments on the deuterated material. [? ] . Fillaux et al's potential is far too broad to be able to describe our 
momentum distribution measurements, pointing up the difficulty of inferring potentials from a vibrational spectrum, 
particularly when the correct assignments are difficult to make. 

The momentum distribution is measured with Deep Inelastic Neutron Scattering(DINS), also called Neutron Comp- 
ton Scattering. DINS is inelastic neutron scattering in the limit of high momentum transfer <f(30 — I00A _1 ). In this 
limit the impulse approximation^ ] can be used to interpret the scattering process which means that neutrons scatter 
from individua protons in the same manner that freely moving particles scatter from each other. Thus the scattering 
function depends only on the proton momentum distribution n(p). In the impulse approximation limit, the neutron 
scattering function is related to the momentum distribution n(p) by the relation 



S(q,u)) = — n(p)S(y - q.p)dp = — J(q,y) (I) 

q J q 

where y = lo — q 2 /2M. S(q,u>) is the Radon transform of n(p). The relation is invertible by a series expansion 
method [? ? ] that represents n(p) in terms of three parameters <7,-, giving the Gaussian widths in three directions, 
and a set of anharmonic parameters, a, n ,i,m, giving the deviations from an anisotropic Gaussian distribution^? ] If 
the & nt i iTn were zero, the potential would be harmonic. 

For temperatures that are well below the lowest possible excitation energy of the proton, n(j>) is determined almost 
entirely by the ground state of the proton, as will be the case with our measurements . The shape of n(p) is 
determined by the Born-Oppcnheimcr potential that localizes the proton. Its second moment gives the kinetic energy 
of the proton. There are in principle corrections to the momentum distribution due to the motion of the oxygen ions, 
but these are insignificant [? ] in harmonic systems as a consequence of the large ratio of the masses of the proton 
and oxygen. We have verified by numerical calculations that they are insignificant also for the anharmonic potentials 
using Path Integral Molecular Dynamics (PIMD). [? ] 

Using a single particle potential to interpret DINS measurements allows us to relate n(p) directly to the proton 
ground state wave function via, 
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We can take the momentum wavefunction to be real if the bond is symmetric, This allows us to calculate the spatial 
wavefunction as well as the BO potential. This is given from the one-particle Schroedinger equation by, 



, , f exp(ip.r) 77— J n(p)dp , 
E - V{r) = i Py F j2mV W_g ^ 

/ exp(ip.r*)-y/ n{p)dp 

For asymmetric bonds, a perturbative approach can be used to extract the symmetric part of the potential surface. 
In this approach, the anharmonic terms in the momentum distribution are related to the anharmonic terms in the 
Taylor expansion of the potential [? ] . 

The DINS experiments were performed on the VESUVIO instrument at ISIS.[? ] Measurements were done on 
single crystal samples of Rb3H(S04)2 at 10K and 70K. The crystal structure for Rb3H(S04)2 is shown in Fig. [TJ 
along with the coordinate system that we will use to describe the potential. The z axis is chosen along the bond, the 
y axis along the line joining the two Rb atoms in the ab plane. Two perpendicular planes of data are collected for 
the crystal sample. One plane is the ab plane where b is the unique axis of the monoclinic crystal. The other plane 
is the ac*, where the c* axis is perpendicular to the ab plane. The ab plane is the plane that contains the Hydrogen 
bond network. The bonds are tilted slightly with respect to it. The unit cell for RbsH(S04)2 has two H bonds at 
60° from each other as shown in Fig. [TJ The momentum distribution of each bond is represented as in Refs.(6,II), 
and the contribution of both bonds added to fit the data. The measurements were done with 30 detectors in a plane, 
distributed with scattering angles from 32° to 68° 

The fitting parameters for the momentum distribution at the two temperatures are given in tables [TJ and [LTJ Many of 
these parameters have a confidence level of 3(7 or higher. There is no phase transition in the protonated Rb3H(S04)2 
for temperatures below 70K. Yet the values of these parameters are accurate enough to reveal some structural changes 
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between 10K and 70K, due simply to the expansion of the crystal. The proton momentum distribution along the 



three axes is shown in Fig. 2(a) Its projection in the yz plane is shown in Fig. 2(b) These figures show no sign of 
proton coherence over two sites, which would show up as an oscillation in the momentum distribution^? ] This can 
be also seen from the BO potential surface plot in the same plane Fig. [3J Along the bond direction, z — axis, the 
potential at both temperatures has a single well with a flat bottom. This is shown more clearly in Fig. 3(a) where 
we show the potential in all three coordinate directions, along with the errors[? Jin the measurement. In general, the 
measurement can only be done out to distances for which the proton has some significant probability of penetrating. 
For our present statistics, this is about 2.5 standard deviations of the spatial wavefunction in the direction of interest. 

Although there is no double well along the bond, it is clear that there is an off axis barrier to motion out of the well 
into the next well. The barrier height is about 350meV, approximately half the value that was inferred from NMR 



measurements on the dcutcratcd material. [? ] at temperatures above room temperature. In Fig. 4(a) we compare 
this potential with the double-well model suggested by Fillaux et a/[? ] to interpret the INS measurements for the 
vibrational spectrum of Rb 3 H(S04)2 at 20A". This model is supposed to produce energy levels that can be assigned 
to the peaks in the INS spectrum. Calculating the mom entum distribution for Fillaux et al's model and comparing 
it with our measurement along the bond direction, Fig. |4(b)[ shows that the inferred potential is in fact far from 
anything that could reproduce our measured momentum distribution, as is also obvious by direct comparison of the 
measured and inferred potentials. Ignoring the motion of the heavy ions, and calculating the excitation energies 
from the measured potential, we find that the first three excitations are at 125 mev, 147 mev, and 155 mev above 
the ground state at 229 mev, with the zero of energy the zero of the potential. There arc only two strong peaks 
identifiable in the INS spectrum, at 142 mev and 192 mev. Fillaux et al assigned these to the transverse modes, and 
identified a band at around 50 mev as the result of the 0-2 transition along the bond interacting with the the SO4 
groups, and a mode at about 5mev as the tunnel splitting. We think in fact that the mode at 192 mev is actually 
a composite of the (coincidentally) nearly degenerate vibration along the bond and the transverse vibration in the x 
direction, which we see at 155 mev and 147 mev respectively. That these should be close in frequency can be seen 
directly from fig 2a and 3a. The lack of accuracy of our predicted frequencies is perhaps due to the size of our error 
bars, which become larger in the region of the excited state energies, where the ground state wavefunction is small. 
Since the excited states probe the potential at distances greater than the ground state does, a fitting procedure that 
combines known spectra with the potential measurements has the possibility of extending the range and accuracy of 
the potential measurements, but we have not attempted that here. 

This single well potential of the short H bond in the protonated Rb3H(S04)2 is the reason behind the absence of 
the low T antiferroelectric transition that exists in the deuterated systems. It length is 2.4A by our measurements 
and 2.483 at room temperature. [? ] The deuterated H bond, on the other hand, 2.518A at 25K [? ], is significantly 
longer than the protonated one, and seems to be, in contrast to the situation in protonated Rb3H(S04)2 , of the 
double minimum type. Rb3H(S04)2 contracts upon cooling which reduces the H bond length. This suggests that 
the double minima potential disappears upon cooling in the protonated system well before the temperature for the 
onset of the AF transition in the deuterated material. This is consistent with ab-initio calculations of the structure 
that take into account the lowering of the energy by derealization of the proton. [? ], and experimental results on the 
related material K3H(S04)2-[? ] 

Double Morse potentials (DMP)have been used to describe O-H-0 bonds phenomenologically in a variety of physical 
systems. [??????] It is of interest to see how close they actually are to the measured potential. We show 
in Fig. 3(a) the best fitting DMP to our data in the yz plane. The fit includes our fit of the 0-0 separation. It can 
be seen that the fit is semi-quantitatively correct. We also tried fitting our momentum distribution measurements 
with a single Morse potential. This would be appropriate if the symmetric crystal structure were only an average, 
and individual bonds were actually off center, covalently bonded to one or the other of the oxygens. That fit doesn't 
work. When we compare the measured momentum distribution with that of a single Morse potential of several 
known models, we see that these models give Gaussian distributions with widths that are significantly bigger than our 
harmonic parameters along the bond direction. For instance, the models by Holzapfel, Matsushita et al, Mashiyama, 
and Borgis et al have momentum distributions of widths 8.2A" 1 , 6.6A -1 , 7.5A -1 and 7.15A -1 respectively compared 
to 4.6A -1 for the harmonic parameter along the bond direction in our measurements. The fact that a symmetric 
DMP fits the measured potential very well and a single Morse potential localizes the proton far too much is some 
indication that the bond is symmetric. The best fit parameters are collected in Tabic Mil and compared with those 
of Holzapfel[? ] and Matsushita and Matsubara [? ]. The parameters were chosen by Matsushita and Matsubara 
to reproduce, in a semi-empirical way, the structural properties of many OH — O complexes. The parameters in 
both works represent the average behavior of wide range of H-bonds. They don't necessarily represent an isolated 
H-bond, as in our measurements. A recent study gives very different parameters for a specific system, KDP[? ], table 
IIII1 which suggests that there is not an accepted set of parameters. However, although the parameters appear quite 
different between our measurement and the parameters of Matsushita and Matsubara, for instance, the potentials, 
over the region in which the proton has significant probability of being found, are very similar, up to an additive 



4 



constant. We show in Fig 4a the comparison of the potential and momentum distribution obtained by Matsushita 
and Matsubara when the 0-0 separation is set to 2AA with no other changes. 

Changing the 0-0 separation in the DMP potential allows for a transition from single to a double well potential 
as the 0-0 separation increases beyond a critical value. The measured critical value for a similar system, K 3 H(S04)2 
is 2AlA [? ] compared to 2.51A in our DMP fit. Comparing our measurements with Ab-initio Calculations made 
for isolated Hydrogen bonds in Aquonium Pcrchlorate [? ], we see that a bond length of about 2.4A is very well into 
the single well regime. The Path Integral Car-Parincllo calculations made by Benoit et al.[? ] give a critical length 
of about 2.36A. According to this calculation a bond of a length of 2AA has a very low barrier, several me Vs. Such 
barriers wouldn't be noticeable in our measurements within the experimental error. This is true as well of the small 
barriers and asymmetry observed in the NMR measurements of Mikacs et al[? ] The single well nature of the bond 
is therefore consistent with all other available evidence. 

We have shown here the first measurements of a 3-D Born Oppcnhcimcr potential in any system. The mea- 
surements can be done on any symmetric hydrogen bond for which the interaction between the protons in dif- 
ferent bonds is negligible. The errors in our measurements are due primarily to counting time limitations on 
VESUVIO, the one existing machine capable of doing NCS experiments. Its count rate could be improved by 
a factor of 10 by the addition of more detectors, and another factor of ten could be achieved by the construc- 
tion of a similar machine on SNS or another newer source. We conclude, therefore, that routine measurements of 
Born-Oppenheimer potentials in hydrogen bonded systems to within ±2meV are feasible with existing technology. 

Acknowledgements: 

The work of G. Reitcr and D. Homouz was supported by DOE Grant DE-FG02-03ER46078. The authors would 
like to thank Phil Platzman for his comments and a careful reading of the manuscript, and Sossina Haile for useful 
conversations . 



■5 



TABLE I: The harmonic scale factors for the momentum distribution measured at 10K and 70K. 

Harmonic Coefficients 

10K 70PT~ ° 
1 <Ji ("A" 1 ) a t ("A" 1 ) 
z 4.60 4.70 
x 4.35 4.36 
y 3/T3 3.65 



TABLE II: The anharmonic fitting coefficients and their variances for the momentum distribution measured at I0K and 70K. 

Anharmonic Coefficients 
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-8.09E-02 


1.0E-02 
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1.0E-02 
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1.49E-03 


9.0E-04 


2.38E-03 


1.01E-03 


10 8 8 


-3.92E-03 


3.0E-03 


-9.48E-03 


2.0E-03 



TABLE III: Comparison of the measured BO potential, fit with a Double Morse potential(DMP) of the form Udmij) = 
U(r — zo) + U(r + zo), where the single morse potentials U(r) = D{exp[— 2a(r — ro)] — 2exp[— a(r — ro)]}, with several other 
DMP models. [Matsushita and Matsubara(MM), Holzapfel(HZ), Mashiyama(MA)] 





This work 


MM 


HZ 


MA 


zo(°A) 
D(meV) 
aCA- 1 ) 
ro(°A) 


1.2±0.07 
2920±205 
1.96±0.01 
0.90±0.06 


2150 
2.89 
0.95 


5337 

2.8±.2 

0.956 


2200 

3.8 

1.00 
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FIG. 1: (a) The A2/a crystal structure for RbzH{S04)'2 and the orientation of our coordinate system. Oxygen(blue), hydro- 
gen(small white), rubidium(olive) , sulphur(red). The gold pyramids are symbols for the SO4 groups. The z — axis represents 
the H-bond direction, the y — axis is perpendicular to the bond in the ab plane, which is the plane with significant high 
temperature proton conductivity. 
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(a) (b) 

FIG. 2: Proton momentum distribution for Rb2,H(SOA)2 at 1QK: (a) along the x — axis (Black), y — axis (Red) and z — axis 
(Green), the rms errors are shown in dashed lines, (b) in the yz plane. 




(a) (b) 

FIG. 3: BO potential for protons in Rb3H(S04)2 at 1QK: (a) along the x — axis (Black), y — axis (Red) and z — axis (Green), 
the curves are shifted by 100 meV along the vertical axis for clarity. The errors are shown in dashed lines, The fit to a double 
Morse potential with parameters shown in the first column of Table III are shown in circles, (b) The potential energy surface 
in the yz plane. 
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FIG. 4: Comparing NCS measurements with the double-well potential model of Fillaux and a Morse Potential Model. Fillaux's 
model is shown in black. The NCS measurements at 10K and 70K are shown in red and green respectively. Matsushita and 
Matsubara's model is shown in blue, (b) The momentum distributions associated with the potentials in (a). The dashed line 
is the fit to a single Morse potential 



